Pseudo-rip: Cosmological models intermediate between the 
cosmological constant and the little rip 
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If we assume that the cosmic energy density will remain constant or strictly increase in the 
future, then the possible fates for the universe can be divided into four categories based on the time 
asymptotics of the Hubble parameter H(t): the cosmological constant, for which H(t) = constant, 
the big rip, for which H(t) — Y 00 at finite time, the little rip, for which H(t) — > 00 as time goes to 
infinity, and the pseudo-rip, for which H(t) —> constant as time goes to infinity. Here we examine 
the last of these possibilities in more detail. We provide models that exemplify the pseudo-rip, which 
is an intermediate case between the cosmological constant and the little rip. Structure disintegration 
in the pseudo-rip depends on the model parameters. We show that pseudo-rip models for which the 
density and Hubble parameter increase monotonically can produce an inertial force which does not 
increase monotonically, but instead peaks at a particular future time and then decreases. 



Introduction — Time is an elusive concept of great in- 
terest to all physicists since it underlies all dynamical sys- 
tems. Although we can measure time on everyday scales 
with exquisite precision, the deeper meaning of time may 
be addressed only by study of the origin and fate of the 
universe. 

Current observations [lpfll strongly suggest that the 
universe is dominated by a negative-pressure component, 
dubbed dark energy. This component can be character- 
ized by an equation of state parameter w, which is simply 
the ratio of the pressure to the density: w = p/p; for ex- 
ample, a cosmological constant corresponds to w = — 1. 
While it is often assumed that w > — 1 in accordance with 
the weak energy condition, it has long been known [B| 
that the observations are consistent with w < — 1, which 
corresponds to a dark energy density that increases with 
time t and scale factor a. If the density increases mono- 
tonically in the future, then the universe can undergo a 
future singularity, called the "big rip," for which p — > 00 
and a — > 00 at a finite time. Shortly before this singular- 
ity is reached, bound structures are disintegrated by the 
expansion Q. 

Note, however, that a dark energy component with a 
monotonically increasing density that is unbounded from 
above does not lead inevitably to a future singularity, al- 
though it does ultimately lead to the dissolution of all 
bound structures. Such models, dubbed "little rip" mod- 
els, were first examined in detail by Frampton et al. 0]> 
who derived the boundary between big rip and little rip 
models in terms of p(a). Properties of little rip models 
were further investigated in Ref. [1]. 

Here we investigate a different set of models, in which 
the density of the dark energy increases monotonically 
with scale factor but is bounded from above by some 
limiting density, poo. Such models can still lead to a dis- 
solution of bound structures for a sufficiently strong iner- 
tial force (which we will define), so we dub these models 
"pseudo-rip" models. This allows us, in the context of 



monotonically increasing p, to distinguish three distinct 
cosmic futures: the big rip, little rip, and pseudo-rip. 
(Note that models for which p(a) is not monotonic are 
physically less plausible, and it is more difficult to make 
any sort of general statement about such models) . While 
the dissolution of bound structures is inevitable in the 
big rip and little rip scenarios, it may or may not occur 
in a pseudo-rip, depending on the model parameters. 

In the next section, we provide the definition of the 
pseudo-rip and examine the conditions necessary to dis- 
solve bound structures. We then examine two specific 
functional forms of the dark energy density and discuss 
scalar field realizations of pseudo-rip models. We assume 
a flat FLRW metric and c = 1 throughout. 

Definition of Pseudo-Rip — Before defining the pseudo- 
rip, we suggest a compact way to classify all rips. We find 
it useful to classify ripping behaviour by means of the 
Hubble parameter H(t). By ripping behavior, we mean 
any future evolution which can lead to disintegration of 
structure in the form of bound systems by virtue of the 
strong inertial force due to dark energy. 

Given the Hubble parameter H(t) for t > to, where 
to is the present time, the density p(t) and pressure p(t) 
are: 
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The big rip is defined by 
H(t) — > +00, 



2H(t)+3H(t)< 



(1) 
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In a big rip, all bound structures dissociate in a finite time 
in the future, and space-time "rips apart" at a finite time 
in the future, i.e., the scale factor of the FLRW metric 



goes to infinity at t = t r i P [6| 
H{t) — > +00. 



. The little rip is denned by 
t -> +00. (4) 



The little rip dissociates all bound structures, but the 
strength of the dark energy is not enough to rip apart 
space-time as there is no finite-time singularity. 

An excellent fit to all cosmological data is provided by 
the ACDM model which, in present parlance, is a "no- 
rip" model defined by 



H(t) = H(t ). 



(5) 



There remains just one additional possibility for mono- 
tonically increasing H(t), namely 



H(t) 



Hoo < CO, t — > +00, 



(6) 



where Hoo is a constant. Equation ^ defines the pseudo- 
rip, the subject of the present article. A pseudo-rip disso- 
ciates bound structures that are held together by a bind- 
ing force at or below a particular threshold that depends 
on the inertial force in the model. Eqs.®-© clearly ex- 
haust all possibilities for a monotonically increasing H(t), 
i.e., a monotonically increasing dark energy density p(a). 

The equations for monotonically increasing H(t) are 
the same for p(t) mutatis mutandis. Pressure p(t) — > 
—p(t) as t — > 00, provided H(t) — > on the R.H.S. of 
Eq.([2|), which is the case for a pseudo-rip model. 

Our division of the future evolution into the categories 
of big rip, little rip, and pseudo-rip represents a differ- 
ent set of models than those examined in Ref. [s| , which 
provided a classification scheme for future singularities. 
Our scheme represents a classification of all models with 
monotonically increasing dark energy density, for which 
the scale factor a goes smoothly to infinity, at either finite 
or infinite time, and for which there are no singularities 
in the derivatives of H unless H itself is singular. In our 
scheme, the big rip encompasses the type I singularity 
of Ref. [§], while the type II, III, and IV singularities 
lie outside the types of models considered in our clas- 
sification scheme. The little rip and pseudo-rip models 
are by definition non-singular, so they fall outside of the 
purview of Ref. . 

The inertial force Fi nert on a mass m as seen by a 
gravitational source separated by a comoving distance I 
is given by [j| 



F 



= ml(H(t) + H(t) 2 ) 

= -ml^-(p(a) + 3p(a)) 

= ml — — (2p(a) + p'(a)a). 



(7) 



For simplicity, we set the scale factor at the present time 
a = a(to) = 1- A bound structure dissociates when 
the inertial force, dominated by dark energy, grows in 
the future to equal the force holding together the bound 



structure in question. For a pseudo-rip, Fi nert is asymp- 
totically finite. 

However, if the bound structure is massive enough to 
significantly affect the local space-time metric, it is not 
accurate to express Fi nert in terms of the FLRW metric. 
A more accurate method and local metric is employed 
in |l0|, and we use their method to calculate the disin- 
tegration times for the Milky Way and the Earth-Sun 
system. 

We analyze two psuedo-rip parameterizations for dark 
energy density, models 1 and 2, each as a function of the 
scale factor a(t) with other parameters. 

Model One — Model 1 is defined by 



Pi(a,B,f,s) 
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where po is the present value of the dark energy density. 
Note that pi is normalized to be po at a = 1, which we 
define to represent the present. Then p\(a,B, /, s) is a 
function of the scale parameter a(t) and of three other 
parameters B, f, and s. It is most convenient to fix / 
and s, which mostly control the strength of the rip, and 
to keep B as a free paramater for fitting the supernova 
data. 

We fix / and s to specify how powerful the psuedo- 
rip should be. The remaining free parameter is chosen 
to make a best fit to the latest sup ernova data from the 
Supernova Cosmology Project ll( with a reduced x 2 of 
-0.98. 

As examples of bound states we consider the Milky 
Way (MW), the Earth-Sun system (ES), the hydrogen 
atom (H atom), and the proton. The first two, MW and 
ES, are gravitationally bound, while for the H atom and 
proton, we must carefully consider the electromagnetic 
and strong color forces respectively. In all cases the dark 
energy density increases monotonically from po at the 
present time to an asymptotic value. Depending on the 
parameters B, /, and s, the inertial force can successively 
disintegrate the MW, the ES, the H atom and the proton. 

In different cases, some of these bounds systems will 
be disintegrated and not others. If none of them are 
disintegrated, we shall refer to such a pseudo-rip as a 
"failed rip." 

In Fig. (pj we show five examples of the scaled Fi nert 
for pi(a), which include the matter and radiation con- 
tributions. Because of these contributions, the curves go 
to negative infinity as a goes to 0, and the ^-intercepts 
are values of a < 1 at which dark energy domination 
begins. Going from bottom to top, the curves represent 
respectively a failed rip; a pseudo-rip which disintegrates 
only the MW; a pseudo-rip which breaks apart the MW 
and ES; a pseudo-rip which destroys the MW, ES, and 
H atom; and the highest curve is for a pseudo-rip which 
succeeds in ripping apart all four of the MW, ES, H atom, 
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FIG. 1: Plotted from the innermost to the outermost curve is 
the scaled Finert for pi (a, 0.29, 10, 48), pi (a, 0.0108, 0.003, 48), 
pi (a, 0.01078, 2 x 10~ 7 ,48), pi (a, 0.01078, 10" 23 , 48), and 
pi (a, 0.0108, 10~ 92 , 48) respectively, for Model 1, given by Eq. 
(JSJ) . Each curve was fitted to supernova data with the ex- 
tra constraint that B > 0. The values necessary for struc- 
tural disintegration are indicated. From the innermost curve 
to the outermost: failed rip; tmw — to = 9.2 x 10 4 Gyr; 
tss - to = 1.6 x 10 9 Gyr; t H atom - to = 1.3 x 10 33 Gyr; 
tproton ~ to — 3.8 x 10 117 Gyr. (The color plots are in the 
online version of the paper.) 



and proton. Note that, unlike a little rip, the asymptotic 
value of piifl) as a — > oo is finite. 

For a given pseudo-rip model with monotonically in- 
creasing dark energy density Pde{p), the structures with 
bigger binding forces disintegrate after those with smaller 
binding forces. But a particular pi can be constructed 
such that it leads to, for example, the disintegration of 
the proton before another p\ disintegrates the Milky Way. 

Notice that the more violent the pseudo-rip is required 
to be, the more extremely small the / parameter is. One 
may counteract this fine tuning of / by, for example, 
introducing new factors into the model that help p\ grow 
faster while still leaving it asymptotically finite. One 
such factor could be (-2±2- ±±M) W . Such a factor is 1 

\a+2q l+q ' 

The 



i+2q l+q ' 

when a = 1, so pi is still po at the present time, 
new parameters q and w can avoid the fine tuning 
Model 2— Model 2 is defined by 



p 2 (a, A, n, m) 



po^-(tan _1 (a — n) — tan _1 (l — n) + 1)' 



(9) 



Like p±, p2 is normalized to be po at a = 1. Then 
P2{a, A, n, m) is a function of the scale factor a(t) and of 
three other parameters A, n, and m. It is most conve- 
nient to fix n and m, which mostly control the strength 
of the rip, and to keep A as a free paramater for fitting 
the supernova data. So we fix n and m to specify how 
powerful the psuedo-rip should be, and the remaining 



free parameter is chosen to make a best fit to the latest 
supernova data from the Supernova Cosmology Project 
11] with a reduced % 2 of ~ 0.98. 

In Fig. ((2|), we plot five examples of the scaled Fi nert 
for P2(a), which include the contributions from matter 
and radiation. Just as in Fig. (fTJ) , the contributions cause 
the curves to approach negative infinity as a goes to 0, 
and the x-intercepts are values of a < 1 at which dark 
energy domination begins. The bottom curve exhibits a 
failed rip, while the others show pseudo-rips of various 
strengths. As was mentioned for Model 1, we see from 
Fig. @ that the disintegration time for the proton for a 
particular p^. can be sooner than the disintegration time 
for the Milky Way for another pi- 

Note in Fig. ([2]) that each Fi nert has a local maxi- 
mum. Because of this bump, p2(a — »• oo) is less than 
the maximum value of p2- So it is possible, for particu- 
lar parameterizations, for Fi nert to reach the level of the 
binding force of a bound structure and then decrease, al- 
lowing the structure to possibly come back together. In 
principle, a pseudo-rip model can have Fi nert with an ar- 
bitrary number of local maxima that give structures the 
chance to dissociate and reform multiple times. All this 
can be achieved using a functional form for dark energy 
density that is monotonically increasing. However, all 
the examples shown in Fig. ([2]) have their asymptotic 
values higher than the values necessary to rip apart the 
structures mentioned in the plot. 
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FIG. 2: Plotted from the innermost to the outermost 
curve is the scaled Fi nert for P2(a, 17481.3, 200, 0.5), 
p 2 (a, 3601.31, 180,2), p 2 (a, 571.1, 160, 10), 

p 2 (a, 171.045, 130,22), and p 2 (a, 55.45, 100, 40) respec- 
tively, for Model 2, given by Eq. @. Each curve was 
fitted to supernova data. The values necessary for structural 
disintegration are indicated. As explained in the text, disin- 
tegration times from different models should not be directly 
compared. From the innermost curve to the outermost: failed 
rip; t M w — to = 45 Gyr; t E s - to — 40 Gyr; t H at om — 1 = 39 
Gyr; t pro ton — to — 38 Gyr. (The color plots are in the online 
version of the paper.) 
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Scalar Field Realizations — One possible realization 
of pseudo-rip models is a minimally-coupled phantom 
model, i.e., one that involves a scalar field with a neg- 
ative kinetic term. The equation of motion for such a 
field is 



4> + 3H<p- V'(4>) = 0, 



(10) 



where the dot is a time derivative, and the prime denotes 
the derivative with respect to <f>. A field obeying this 
equation of motion rolls "uphill" in the potential. 

It is clear that a sufficient condition for a pseudo-rip 
is that V(4>) — > Vq (where Vo is a constant) as <j> — >• oo, 
and in Ref. Q it was shown that this is also a neces- 
sary condition for a monotonic potential. In this case, 
we simply have p^ = Vo- If the potential is not mono- 
tonic, the density of the scalar field can also approach a 
constant asymptotically if the field gets trapped in a local 
maximum with V — Vq. Phantom fields with bounded 
potentials have been discussed previously in Ref. [12 1. 

Note, however, that our discussion of pseudo-rips is 
much more general than the specific example provided by 
phantom field models. Phantom fields represent only a 
single possible realization of this much more general class 
of models for the asymptotic expansion of the universe. 

Discussion — We have described merely two illustrative 
models of the pseudo-rip. Obviously, there is an infinite 
number of possibilities. 

A failed rip will disintegrate nothing because the iner- 
tial force will not reach a high enough magnitude. Model 
1 exemplifies a pseudo-rip model which can variously dis- 
integrate between one and all four of the chosen systems 
while pi asymptotes to a finite density, which implies a 
finite inertial force, as the scale factor a(t) approaches 
infinity. As can be seen from Fig. (fJJ), these quantities 
plateau to a constant value after the last disintegration 
has taken place. 

We included Model 2 as a particularly interesting ex- 
ample in which the inertial force rises relatively abruptly 
before it plateaus, as shown in Fig. @. The various 
disintegrations take place soon after each other in a cos- 
mological sense. 

For a function that asymptotically approaches a con- 
stant, such as p\{a) or /92(a), a point of inflection al- 
lows the function to have an arbitrarily large slope for a 
portion of the domain and still increase monotonically. 
A point of inflection is present in all the parameteriza- 
tions of pi (a) and p2 (a) and their derivatives in this letter 
and it allows them to fit the supernova data (which re- 
quire the densities to have a very small slope over the 
relevant portion of a) and still reach a high density in 
a relatively short time. However, the relevant quantity 
that determines when a structure dissociates is Fi nert , 
which is proportional to 2p(a) + p'{a)a. This combi- 
nation of terms is responsible for any local maxima in 
Fi ner t, not merely p or p'. An example of a model that 
has an inflection in p' DE {a) but none in pde{&) is given 



by Pde{o) = ct(a — ln[l + e a ~ c ]), where a, C > are con- 
stants. Such a model has a local maximum in Fi nert and 
if we shift by D, where D > C is a constant, neither the 
resulting pde{& + D) nor its derivative has an inflection 
point, but the resulting Fi nert still has a local maximum. 

It is amusing to take examples of Model 2 which are 
more extreme than the ones illustrated. It is possible 
to design a pseudo-rip model such that disintegrations 
happen arbitrarily soon after the present time while still 
maintaining excellent fits to the supernova data. The 
Sun may not rise tomorrow. 

This is a dramatic illustration of the fact that any 
amount of observational data, necessarily restricted to 
the past lightcone and necessarily with non-zero errors, 
cannot predict anything mathematically about the future 
even one hour hence without further assumptions. It is 
also a display of the difference between mathematics and 
physics: the physicist necessarily employs intuition about 
the real world. 

The earliest support for cosmological futures such as 
the big rip, little rip, or pseudo-rip might come from the 
Planck satellite, which is expected to have data around 
September 2012. If the dark energy equation of state 
emerges with w < —1, it will be a shot in the arm for 
such exotic ideas. 

It will also lend new understanding of the nature of 
time and perhaps the beginning and possible cyclicity of 
the universe. 
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